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Several thermodynamic measurements for the cobaltate superconductor, Na^CoCh • 1/H2O, 
have so far provided results inconsistent with each other. In order to solve the discrepancies, we 
microscopically calculate the temperature dependences of specific heat and superfluid density 
for this superconductor. We show that two distinct specific-heat data from Oeschler et al. and 
Jin et al. are reproduced, respectively, for the extended s-wave state and the p-wave state. 
Two different superfluid-density data are also reproduced for each case. These support our 
recent proposal of possible two different pairing states in this material. In addition, we discuss 
the experimentally proposed large residual Sommerfeld coefficient and extremely huge effective 
carrier mass. 
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The Co-oxide superconductor Na^CoC^ -2/H 2 has at- 
tracted great interest since its discovery. 1 Several exper- 
iments have suggested an unconventional superconduc- 
tivity due to a spin-fluctuation mechanism. 2 

A number of thermodynamic measurements have been 
performed to elucidate its superconducting (SC) state. 
However, the reported results on the upper critical 
field (H C 2), specific heat (C(T)) and superfluid density 
(p s (T)) are rather scattered and contain a lot of dis- 
crepancies. As for H C 2, some groups reported values ap- 
parently exceeding the Pauli limit of the weak-coupling 
theory, while the other groups reported values compara- 
ble to, or less than the limit. 2 As for C(T), Jin et al. 
observed a shoulder-shaped structure accompanied by a 
steeper decrease at low temperatures, 3 while Yang et al. 
and Oeschler et al. did not observe such a structure. 4 ' 5 
As for p s {T), Uemura et al. and Kanigel et al. reported a 
linear evolution with decreasing temperature from /xSR 
measurements, 6 ' 7 while Yuan et al. observed a rather 
unusual behavior with negative curvature as shown in 
this paper. 8 In addition, the former two groups proposed 
unusually huge effective carrier mass. 6 ' 7 

On the other hand, several experiments have recently 
revealed that the Co02-layer thickness of this material 
varies sensitively depending on the Na content and the 
amount of H2O molecules, and, furthermore, that the 
magnetic and SC properties are quite sensitive to the 
layer thickness. 2 ' 9-17 Motivated by these findings, we 
previously studied how the variation of Co02-layer thick- 
ness affects the properties of this material. 18 ' 19 In these 
studies, we pointed out the following aspects; [1]: The 
slight layer-thickness variation induces drastic changes in 
the band structure and the Fermi surface (FS) topology. 
In particular, the FS is constructed from double ai s -band 
cylinders around the T-point in the thick-layer systems, 
while in the thin-layer systems, it is constructed from 
one a\ g cylinder around the T-point and six e' g pockets 



near the K-points. The former type of FS is referred to 
as FS1 while the latter is FS2. [2]: The spin fluctuation is 
enhanced at q ^ in the thick-layer systems and at q = 
in the thin-layer systems. 19 The critical enhancement to- 
ward magnetic ordering (MO) is expected in the sys- 
tem with moderate thickness. In addition, T-dependence 
of the spin fluctuation also has strong layer-thickness 
dependence. The discrepancies of the NMR, NQR and 
/iSR results are explained by this layer-thickness depen- 
dence. 19 [3]: Both singlet extended s-wave pairing 20 and 
triplet p-wave pairing are possible in this material de- 
pending on the Co02-layer thickness or the FS geome- 
try. The former is favored in the thick-layer systems with 
FS1, and the latter in the thin-layer systems with FS2 
as shown before. The experimental phase diagram 2 ' 9 ' 10 
containing successive SC, MO and another SC phases is 
explained by this FS deformation. 

In this letter, we discuss the thermodynamic proper- 
ties of the SC states in Na2,Co02 • 2/H2O, and show that 
the scattered experimental results can be explained fairly 
well if we consider the two different SC states with differ- 
ent FS topologies. Certainly, the scattered H c2 data are 
easily understood if we attribute the data exceeding the 
Pauli limit to the triplet p-wave pairing, and the other 
data to the singlet extended s-wave pairing. On the other 
hand, C(T) and p s (T) are quantities sensitive not only 
to the SC gap structure but also to the electronic struc- 
tures such as FS topology and density of states (DOS). 
Thus, in the following, we microscopically calculate these 
quantities (1) for the extended s-wave state on FS1, and 
(2) for the p-wave state on FS2. We show that the two 
different C(T) data reported by Oeschler et al. and Jin 
et al. are reproduced for each pairing state. Two distinct 
behaviors of p s (T) data are also reproduced for each case. 
We also discuss the experimentally proposed large resid- 
ual DOS and extremely huge effective mass of electrons. 
These results support our theoretical proposals about the 
two kinds of FSs in the present cobaltate system. 
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Fig. 1. Band dispersions near the Fermi level. Two types of the 
Fermi surfaces, FS1 (FSla and FSlb) and FS2 (FS2a and FS2b), 
predicted in the previous study are reproduced. 18, 19 



The calculations are performed in the same way as 
ref. 21. We start with the multiorbital Hubbard model: 
H mo = -ffkin. +-ffint.- The first term -f/kin. expresses elec- 
tron hoppings between the Co t2 g orbitals. In this term, 
we consider up to the third nearest neighbor hoppings, 
which are represented by twelve kinds of transfer inte- 
grals ti, tu and a crystal field (CF) parameter A for 
the trigonal CF from O ions. The detailed expression of 
-Hkin. appears in ref. 19. The second term -Hint, is the 
Coulomb-interaction term, which consists of five contri- 
butions as fli n t. = Hu+Hu>+Hj h +Hji + H v . Here, Hu 
and Hjji are the intra- and inter-orbital Coulomb repul- 
sions, respectively, and Hj H and Hj> are the Hund's-rule 
coupling and the pair hopping, respectively. These inter- 
actions are expressed using Kanamori parameters, U, U' , 
Jh and J', which satisfy the relations; U' — U — 2 Jh and 
Jh = J'- The last term Hy = VY], „■ rijra,- expresses the 
Coulomb repulsion between electrons on the adjacent Co 
sites. 

The values of t\, t\2 and A are determined so as to 
reproduce the band structures which predicted the two 
types of FSs, i.e., FS1 and FS2. 18 ' 19 The obtained values 
are listed in Table. I where the energy unit is meV. Fig- 
ure 1 shows band structures with FS1 (FSla and FSlb) 
and those with FS2 (FS2a and FS2b) near the Fermi 
level. The inner a% g cylinder of FSlb (the e' g pocket of 
FS2b) is slightly larger than that of FSla (FS2a). The 
whole Co t2 g band structures are shown in refs. 18 and 
19. Hereafter, the energy unit is eV. 

We derive the effective pairing interactions for both 
spin-singlet and spin-triplet channels (T s (q) and r*(q)) 
by using the random-phase approximation (RPA). Then, 
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Table I. Transfer integrals and crystal-field parameters. The en- 
ergy unit is meV. 




Fig. 2. Momentum dependence of the gap functions for (a) the 
extended s-wave state on FS1, </>„ s (k), (b) the p^-wave state 
on FS2, $T(k), (c) the p H -wave state on FS2, </>£ a (k), and (d) 
the linear combination of the p^-wave and the py-wave states, 

0g(k) = y^WWF+WW- 



we construct the Eliashberg equation. The detailed ex- 
pressions of T s (q), r*(q) and the Eliashberg equation 
are included in refs. 19 and 22. By solving the Eliash- 
berg equation, we obtain the momentum dependence of 
SC gap function </> mn (k). The calculations are numer- 
ically carried out by taking 128x128 k- meshes in the 
first Brillouin zone and 1024 Matsubara frequencies. We 
take 17=0.50, J H =0.05, U=0.05, fc B T=0.01 and the Co 
valence s =+3.4. Using the unitary matrix t/(k) which 
diagonalizes -ffkin., we obtain the gap function with band 
basis ((/> a (k)) from those with orbital basis (</> mn (k)). In 
Fig. 2, we display Q (k) for (a) the extended s-wave state 
on FS1 (C(k)), (b) the p x -wave state on FS2 (^(k)), 
and (c) the p y -wave state on FS2 (^(k)). 

As for the p-w&ve pairing, the p x and p y states are 
degenerate in the triangular lattice. In addition, there are 
three choices of the direction of the d vector (x, y, z), i.e., 
the S z component of the S = 1 Cooper pair. Below T c , a 
linear combination of these six states should be realized, 
namely, z (p x ± ip y ) , xp x ± yp y and xp y ± yp x . All of these 
gaps are represented as <^(k) = y/ ^^(k) 2 + </>a y (k) 2 , 
which is shown in Fig. 2 (d). 

In the following, we assume that the k- and In- 
dependence of the SC gap function below T c is described 
by A kcl (T) = a (k)A(T). Substituting this into the stan- 
dard BCS gap equation, 

1 tanh[ f, k '^ l 
A ka (T) = -- V^ aM , a -^llAva<(T), (1) 

k',a' ' k 



w ith Vka,k'a' = -V0 o (k)0*/(k') and Eu a 

VZL + \ A **( T )\ 2 > we obtai ^ 

i = -^Ei*.wi 2tani[fel 



(2) 



Considering that A(T c )=0, we calculate the tempera- 
ture evolution of A(T) below T c by solving eq. (2). 
Along with T c of the actual material, we take fcBT c =0.45 
meV (T c ~ 4.5 K) throughout our calculations. The en- 
tropy S(T), C(T) and p s (T) are calculated by the for- 
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Fig. 4. Calculated superfluid density p s (T)/p s (0.045T c ) and ex- 
perimental u s (T) data from Uemura et al.® and Kanigel et al.J 
and A(T) — 2 data from Yuan et ai. 8 Black dashed lines are 
crude extrapolations of the data to T — » 0. Inset: Calculated 
p s (T)/p s (0.045T c ) (o- s (T)/cr s (0.045T c )) for the p-wave state on 
FS2b. The contribution from the a\ g FS and that from the e' g 
pockets are shown separately. 



Fig. 3. Calculated specific heat C(T)/^ n T for (a) the extended 
s-wave state on FSla and (b) the p-wave state on FS2a. The ex- 
perimental data from refs. 5 and 3 are also depicted respectively 
in each figure. A crude extrapolation gives a superficialy large 
7^ as shown in (a) by the dashed line, (c) Momentum depen- 
dence of the gap functions of both pairing states, </>^ s (k) (red 
line) and </>2(k) (blue line), together with the band dispersions 
near the Fermi level for both FSla and FS2a cases (red and blue 
dashed lines, respectively). 
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Figures 3 (a) and 3(b) show the calculated Sommer- 
feld coefficient r y(T)—C(T)/T normalized by the normal- 
state 7 ( 7n ) for (a) the extended s-wave state on FSla 
and (b) the p-wave state on FS2a, together with the ex- 
perimental data from Oeschler et al. (ref. 5) and Jin et 
al. (ref. 3), respectively. Both figures show good coin- 
cidece between the experiment and the calculation. In 
the following, we discuss the results in detail. 
[1]: In Fig. 3(a), we show the contributions from the two 
parts of FS1. Apparently, the inner FS dominantly con- 
tributes to the jump at T c (A7) for the FS1 case. This is 
because the extended s-wave state on FS1 has a large gap 
amplitude on the inner FS. This can be seen in Fig. 3(c), 
in which ^>„ s (k) is depicted together with the band dis- 
persion for the FSla case near the Fermi level. This fig- 
ure also shows that the gap on the outer FS of FS1 is 
markedly small. Similarly, for the p-wave state on FS2, 
the cylindrical a\ g FS dominantly contributes to A7 as 
shown in Fig. 3(b). This is due to a large gap ampli- 
tude on the ai g FS as shown in Fig. 3(c). In the present 
p-wave state, the gaps on the e' g pockets are markedly 
small. Note that this is in sharp contrast to the previous 
proposals of the p-wave or /-wave state with large gap 
amplitudes on the e' g pockets. 22-25 

[2]: For the p-wave state on FS2, the gap almost van- 
ishes at certain points on the e' pockets, which form 



line nodes. Consequnce of such a line-nodal gap is seen 
in the T 2 -dependence of C(T) ( 7 (T)=C/T cx T) at low 
temperatures in Fig. 3(b). 

[3]: Oeschler et al. reported H c2 (T -> 0) of -8 T, which 
is comparable to, and does not exceeds the Pauli limit of 
~8.3 T for T c ^4.5 K. 5 This is in agreement with the spin- 
singlet pairing expected in their sample. On the other 
hand, Jin et al. reported H c2 (T -> 0) of 17.1 T, which 
apparently exceeds the Pauli limit. 3 This is consistent 
with the expected spin-triplet pairing in their sample. 
[4]: The C(T) data of Jin et al. shows a deviation 
from linearity with steeper decrease below 0.8 K 
(T/T c =0.17 and T c =4.7 K). Our result well reproduces 
this behavior as shown in Fig. 3 (b). In contrast, Oeschler 
et al. did not observe such a decrease down to 0.85 K 
(T/T c =0.19 and T c =4.52 K). However, as shown in Fig. 3 
(a), if the temperature is further decreased, the steep de- 
crease should be observed also for the extended s-wave 
case. 

[5]: Oeschler et al. claimed a large residual 7 (7 re s) of 
6.67 mJ mol _1 K -2 at T — > by extrapolation. The esti- 
mated SC volume fraction (j n — 7 rC s)/7ri is small to be 
59 % where 7n =16.1 mJ mol -1 K -2 . In fact, a number 
of measurements have reported similar or larger 7 rcs so 
far. 2 Our results suggest two reasons why the experimen- 
tally estimated 7 rcs is so large. One is possible existence 
of non-SC sections of the FS. Both extended s-wave and 
p-wave SC gaps have markedly small amplitude on the 
outer FS of FS1 and on the e' g pockets of FS2, respec- 
tively. These tiny gaps should be easily destroyed by im- 
purities and/or defects. The resulting non-SC sections of 
the FS should generate large 7 ros . In fact, this possibility 
was also proposed by Takada et al. 2 A large residual DOS 
reported from NQR 59 (l/TiT) measurements 26 may also 
be attributed to such non-SC FSs. Another reason is the 
crude extrapolation which overlooks the steep decrease 
at very low temperatures. This should cause the overesti- 
mate of 7 res as indicated by the dashed line in Fig. 3(a). 
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Next, we show the calculated p s (T) in Fig. 4 for (a) 
the extended s-wave state on FSlb, and (b) the p-wave 
state on FS2b, which are normalized by the value at 
T=0.045T C . In addition, the results for the uniform BCS 
gap are shown by dashed lines for comparison. Both the 
extended s-wave and the p-wave cases exhibit unusual 
behaviors, which are completely different from the BCS 
case. Below T c , both two start to increase rather gradu- 
ally, and even at low temperatures, both exhibit increas- 
ing behaviors in contrast to the saturating behaviors of 
the BCS cases. We compare our calculations with the 
experimental data of pSR rate a s (T) reported by Ue- 
mura et al. (ref. 6) and Kanigel et al. (ref. 7), and the 
magnetic-penetration-depth (A(T) -2 ) data measured by 
Yuan et al. utilizing a tunnel diode based self- inductive 
technique (ref. 8). Here, <r s (T) and A(T) are related to 
p s {T) as a s (T) oc A(T)~ 2 = 4irp s (T)e 2 /m*c 2 where m* 
is the effective mass of the carriers. We also depict their 
data with proper normalization in Fig. 4. 

For the extended s-wave case, the calculated p s (T) 
(blue line) is convex downward, and shows quite good 
coincidence with the experimental data from Yuan et al. 
The H C 2 value measured in this experiment is less than 
the Pauli limit, which is consistent with the present spin- 
singlet pairing. For the p-wave case, the calculated p s (T) 
(red line) shows an steeper increase starting from a lower 
temperature than the extended s-wave case and shows 
good agreement with the expreimental data from Kanigel 
et al. and Uemura et al. However, the data from Uc- 
mura et al. slightly deviates from the calculation, which 
shows no steep increase at low temperatures. In fact, 
the steep increase is originated from the weakly gapped 
e' g pockets (See the inset. Contributions from each part 
of FS2 are depicted separately). Such small gaps would 
be easily destroyed by impurities and/or defects, result- 
ing in the residual normal-state carrier density p n . This 
should cause the suppression of the increasing behavior 
ofp s (T). 27 

Both Kanigel et al. and Uemura et al. estimated <r s (T) 
at T — > by extrapolation. They assumed that thus ob- 
tained cr° xt (0) corresponds to a s (0) for the clean limit 
(tTg lean (0)) where all the normal-state carriers above T c 
contribute to the superconductivity at T —0 (namely, 
p s (T = 0)=p n (T > T c )). Using this of oan (0) value, the 
effective carrier mass m* was estimated to be about 75 
times the bare electron mass m e by Kanigel et al. and 
about 100 times by Uemura et al. These values are un- 
usualy huge and comparable to the effective mass in the 
heavy-electron superconductors. However, the extrapo- 
lation leads to underestimate of of ean (0) when we over- 
look the subsequent steep increase at low temperatures 
or the possible large residual p n (0) arising from the non- 
SC e' g pockets due to the sample inhomogeneity. As a 
result, their extrapolations give only 13-20 % of the ac- 
tual value of af can (0). (Note that in Fig. 4, the crude 
extrapolations give cr° xt (0)/CT s (0.045T c ) ^0.6 for Kanigel 
et al. and ^0.4 for Uemura et al. wheareas the calcu- 
lation gives of ean (0)/o- s (0.045T c )=3.06.) These under- 
estimates should cause approximately 5 and 7.5 times 
overestimates of m*, respectively. Thus, the estimated 
effective mass should be m* ~15 m e . This value is rea- 



sonable, and is in good agreement with that estimated 
from the angle-resolved photoemission spectroscopy. 28 

To summarize, we have studied the T-dependence of 
C(T) and p s (T) in the two SC states of Na^CoOa -yH 2 
predicted in our previous calculations; the extended s- 
wave state on FS1 and the p-wave state on FS2. We 
have shown that the two different C(T) data as well 
as the two distinct p s (T) data are reproduced for the 
two cases. This advocates our proposal of possible two 
different SC states in this material depending on the 
Co02-layer thickness or the FS geometry. In addition, 
there are weak SC sections of the FS for both cases; the 
outer FS of FS1 and the e' g pockets of FS2, respectively. 
These cause sharp decrease (increase) of C(T) (p s (T)) 
at very low temperatures, which makes the extrapola- 
tion of the data to T — > difficult, and generates large 
residual 7 and p n when these small gaps are destroyed by 
impurities and/or defects. The reported unusually large 
7 ros and the extremely huge m* are understandable from 
these facts. We lastly note that it is quite possible that 
the two superconducting states can coexist in the same 
sample, especially in large single crystals, due to the un- 
uniformity of the Co02-layer thickness or the ununiform 
water distribution, as suggested from a recent impurity- 
effect measurement. 29 
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